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ABSTRACT 

Context. The calibration of binary systems with accurately known masses and/or radii provides powerful tools to test stellar structure and 
evolution theory and to determine the age and helium content of stars. We study the eclipsing double-lined spectroscopic binary system RS 
Cha, for which we have accurate observations of the parameters of both stars (masses, radii, luminosities, effective temperatures and metallicity). 
Aims. We have calculated several sets of stellar models for the components of the RS Cha system, with the aim of reproducing simultaneously 
the available observational constraints and to estimate the age and initial helium abundance of the system. 

Methods. Using the CESAM stellar evolution code, we model both components starting from the initial mass and metallicity and adjusting the 
input parameters and physics in order to satisfy the observational constraints. 

Results. We find that the observations cannot be reproduced if we assume that the abundance ratios are solar but they are satisfied if carbon 
and nitrogen are depleted in the RS Cha system with respect to the Sun. This is in accordance with the abundances observed in other young 
stars. The RS Cha system is in an evolutionary stage at the end of the PMS phase where models are not strongly sensitive to various physical 
uncertainties. However we show that the oscillations of these two stars, which have been detected, would be able to discriminate between 
different options in the physical description of this evolutionary phase. 
Conclusions. 
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1. Introduction 

Stellar calibration that-is determining masses, radii, ages of 
stars using theoretical stellar evolutionary tracks is a power- 
ful tool. However, its results strongly depend upon the valid- 
ity of the adopted stellar models. The physical description of 
these models can be validated with a small set of stars whose 
masses, radii, luminosities, effective temperatures and metal- 
licities are accurately known. An eclipsing binary system is 
therefore an excellent candidate for such a test as the orbital 
information provides masses and radii and the assumed com- 
mon origin of both components (implying same age and chem- 
ical composition) brings out additional severe constraints to the 
modelling. Validation of main sequence (MS) stellar evolution 
models has been quite extensively performed (e.g. Noels et al. 
1991, Morel et al. 2000). On the other hand, only a few valida- 
tion works have sofar concerned the pre-main sequence (PMS) 
phase. Palla & Stahler (2001) tested their PMS theoretical evo- 
lutionary tracks using eight binary systems, assuming a solar 
metallicity in all of them, although this assumption appeared 
later on not valid in all but one case. Marques et al. (2004) mod- 
elled the binary system EK Cep using a ^-minimisation and 
obtained the most reliable set of theoretical stellar parameters 
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which reproduces the observational one. These authors suc- 
ceeded in modelling the PMS secondary component but failed 
to obtain a theoretical radius which reproduces the observed 
one for the MS primary component. 

The A-type system RS Chamaeleontis (RS Cha) is an 
eclipsing double-lined spectroscopic binary system. Mamajek 
et al. (1999) reported X-ray emission from the rj Cha cluster, 
which indicates the PMS status of both components of RS Cha. 
Up to recently, all fundamental parameters of the components 
of this system were known, except the metallicity. However, 
the knowledge of the metallicity is crucial to model the struc- 
ture and the evolution of these PMS stars (see Sect.4.3.2). The 
metallicity of RS Cha system which is used here has been ob- 
tained with spectroscopic data collected at the SAAO (South 
African Astronomical Observatory) (Alecian et al., 2005; here- 
after paperl). We took advantage of these data to redetermine 
the masses and radii of both stars. The resulting physical pa- 
rameters of both components are detailed in Tabled 

We are then able to model these two stars and confront the 
results to observations. This paper presents an advanced study 
of the modelling of the pre-main sequence RS Cha system. The 
next section describes the physical input of our standard evo- 
lution models. Sect. 3 states the observational constraints that 
our RS Cha models must reproduce. As our standard models 
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Tabl e 1. Fundamental pa rameters of RS Ch a. R00: 
iRibas etalJ d2000h. CN80: IClausen & Nordstroml Jl980l> . 
MOO: iMamaiek etalJ d200(A . P stands for primary and S for 
secondary 



P S References 



M/M 


1.89 ±0.01 


1.87 ±0.01 


paperl 


R/Re 


2.15 ±0.06 


2.36 ± 0.06 


paperl 


7W (K) 


7638 ± 76 


7228 ± 72 


R00 


log(L/L ) 


1.15 ±0.09 


1.13 ±0.09 


L = 4nR 2 crT^ 


log(g)(cm.s- 2 ) 


4.05 ± 0.06 


3.96 ± 0.06 


g = MG/R 2 ' 


v sin i (km.s -1 ) 


64 ±6 


70 ±6 


paperl 


P (day) 




1.67 


paperl 


i(°) 


83.4 ± 0.3 


CN80 


[Fe/H] 


0.17 ±0.01 


paperl 



cannot reproduce the whole set of observations, Sect. 4 inves- 
tigates the sensitivity of PMS models to uncertainties in the 
input physics and physical parameters. Sect. 5 presents the final 
calibrated models for the RS Cha system which reproduces all 
the available observations. Sect. 6 contains an analysis of the 
recently discovered pulsations in both components of RS Cha 
and discusses theoretical oscillation modes and periods. Finally 
conclusions are presented in Sect. 7. 

2. Stellar modelling 

We consider both stars as isolated stars since the system is de- 
tached. Unless it is specified, the temporal evolution of the in- 
ternal structure of both stars is carried out at constant mass; 
hence the initial mass is set equal to the observed one. Diffusion 
is negligible over the time duration of a PMS phase, hence we 
assume that the observed surface abundances correspond to the 
initial chemical composition. We stop the evolution when the 
calculated effective temperatures and luminosities correspond 
to the observed ones within their error bars. We stress that in 
the case of a binary system, one additional constraint exists: the 
age of both components must be the same. 

O ur evolutiona ry models are computed with the CESAM 
code jMorelll997l) : they are standard in the sense that no effect 
of magnetic field is included. Effects of rotation and diffusion 
are not included as the time scale of these phenomena are very 
long compared with the time spent by a star on its PMS (20Myr 
for a 2 M ). 

We used the OPA L equation of state jRqgers et alJ fl996) 
an d OPAL opacities Jlglesias & Rogerdll9 96) complemented 
by lAlexander & Ferguson! i 19941) opacities at low tempera- 
tures. 

During its PMS phase, a star can possess either a convec- 
tive envelope or a convective core or both. The temperature 
gradient in convection zones is computed using the classical 
mixing-length theory. The mixing length is defined as I = aHp, 
a being the mixing length parameter and Hp the local pressure 
scale height, Hp = --^p- We set the a value equal to 1.62 as 
obtained for a calibrated solar model calculated with the same 
input physics. We define the overshooting d ov as the length of 
penetration of the convective eddies into the adjacent radia- 
tive zone: d m , = a or H P with a ov the overshooting parameter. 
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Fig. 1. Evolutionary tracks for 1.6-2.2 M standard models in 
a HR diagram. Crosses represent observational error bars in ef- 
fective temperature and luminosity (Table[0 for the primary (P) 
component (left) and for the secondary (S) component (right) 
of RS Cha. Dots indicate the onset of the radiative core, dia- 
monds signal the convective envelope disappearance and trian- 
gles mark the convective core apparition. 

Unless otherwise specified, our models are computed without 
overshooting, i.e. d ov = 0. 

The species entering the nuclear network are: 'H, 3 He, 4 He, 
I2 C, I3 C, 14 N, I5 N, I6 0, n O. We considered 2 H, 7 Li and 7 Be in 
equilibrium and only the most important reactions of PP+CNO 
cycles are taken into account. The nuclear reaction rates are 
taken from NACRE compilation (Angulo et al. 1999). 

The metallicity defined as [Fe/H] = log( ,5/^ ) where Z, 
X are the metal and hydrogen mass fractions respectively, is 
an observational constraint. Its initial value is taken equal to 
the observ ed one ( Table H) and the abundance ratios are the 
solar ones JGrevesse & Noelsll993l) . We take the initial helium 
mass fraction equal to the value from the solar calibration: Yq = 
Y Q = 0.267, and the hydrogen mass fraction is given by Xq = 
i+(z/x) • ^ e atmos phere is calculated assuming the Eddington 
grey law. 

Each evolution is initialized with a homogeneous, fully 
convective model in quasi-static contraction and we define the 
age of a star as the time elapsed since initialisation. Although 
our initial models are not realistic, recent studies have com- 
pared the evolutionary tracks during the PMS phase, consider- 
ing two different initial conditions: the first model being cal- 
culated from the initial cloud contraction and the second one 
being calculated only from the Hayashi track (Cariulo et al. 
2004, F. Palla, private communication). No fundamental differ- 
ences between the two models were found and ages calculated 
in both cases are almost identical. So the initial conditions for 
the formation of the PMS star do not interfere in our study. 

These stellar models will be considered from now on as our 
standard models. 

Fig. H shows PMS evolutionary tracks from the Hayashi 
line to the ZAMS for masses in the range of interest for RS 
Cha. We have indicated by diamonds the locus where the star 
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Fig. 2. Evolutionary tracks of standard models (dash-dotted 
lines) and error boxes in masses and radii for the primary (P) 
and the secondary (S) components. Crosses represent observa- 
tional error bars in effective temperature and luminosity (Table 
^ for the primary (P) component (left) and for the secondary 
(S) component (right) of RS Cha. Lowest to upper tracks are 
for 1.86, 1.87, 1.88, 1.89 and 1.90 M . 

becomes totally radiative and by triangles the locus where a 
convective core appears in the center of the star. At the stage 
where RS Cha is encountered, a star, evolving from the Hayashi 
line towards the ZAMS, has developped a fully radiative enve- 
lope and a small convective core. The earliest nuclear reactions 
of the CNO cycle and the p-p chain begin. The energy released 
by the C l2 (p,y)N n (fi + , v)C 13 (p, y)N 14 chain overheats the in- 
ternal regions of the star. The energy excess builds an excess 
pressure gradient which slows down the contraction rate in the 
center. At the same time, the overheating leads to the apparition 
of a convective core. The central regions begin to expand. The 
luminosity then decreases because not only do these expanding 
internal shells not contribute any more to the energy flux, but 
they blanket the nuclear source, by absorbing a large fraction 
of the nuclear energy produced (see detailed description of the 
PMS evolution in Iben (1 19651) ). For later discussions, it is im- 
portant to note that the primary, as it is slightly more massive, is 
in a more advanced stage of evolution than the secondary one. 
As a result, its luminosity has already decreased unlike that of 
the secondary one (Fig.[2j. 

3. Observational contraints in a HR diagram 

Table \l\ details the physical parameters of both components 
of RS Cha derived from observations. Effective temperatures 
and luminosities as well as masses and radii will then be com- 
pared to the corresponding calculated parameters obtained us- 
ing standard models. We have calculated evolutionary tracks 
for the observed masses of both components as well as for their 
extreme values given by the error bars. Fig|2]shows the end of 
the PMS phase of these tracks for 1.86, 1.87, 1.88, 1.89 and 
1.90 M Q masses. 

For a given mass M t, s , a stellar model is evolved till its 
radius matches the observed radius R i, s . This provides the the- 



oretical luminosity and effective temperature (L, T e ff)caic of the 
model. The same is done for the extreme values of the observed 
masses and radii and this gives rise to error boxes which are 
superimposed onto the tracks in Fig|2]in full lines. In the same 
graph, crosses represent the error bars in luminosity and ef- 
fective temperature for the primary (in left) and secondary (in 
right) components derived from observations (L, T e ff) bs- 

When both crosses are located inside their respective error 
boxes, the corresponding stellar models reproduce the observa- 
tions. As it can see in Fig. |2 this does not occur when stan- 
dard models, built as detailed in the previous section, are used. 
Moreover, the calculated ratio of the primary to the secondary 
luminosities, (Lp/Ls) ca ic, is lower than one, contrary to the ob- 
servations, (Lp/Ls) bs - 1-05 (Table Our standard models 
indeed indicate that the luminosity of the P component has be- 
gun to decrease, unlike the luminosity of the S component. This 
implies that the P component has started nuclear burning earlier 
than the S component (see Sec. 2). 

In order for the numerical models and the observations to 
agree, in particular to invert the (Lp/Ls) ca ic ratio, the error box 
of the primary ought to be shifted towards greater luminosities 
with respect to the secondary box. According to the evolution- 
ary phase of the primary as described in the previous section, 
this means that it is necessary to delay the onset of the CNO 
cycle for the primary. 

The way this can be achieved is carried out in the next sec- 
tion where we investigate the sensitivity of the stellar models 
in this evolutionary phase to various parameters and physical 
inputs. 



4. Sensitivity to physical inputs and parameters 

4.1. Unchanged error boxes in a HR diagram 

Some physical inputs and parameters of stellar model shave no 
effect on the location of the error boxes of RS Cha in a HR 
diagram: 

4.1 .1 . Mixing length and overshoot parameters: 

An increase in the mixing length parameter a moves a PMS 
evolutionary track towards larger effective temperature and to- 
wards slightly larger luminositi es, when the s tar possesses a 
convective envelope (see B ohm-Vitensel {T992)). But, as seen 
in Fig^ the RS Cha stars have no longer a convective enve- 
lope and changing a therefore has no effect on the error boxes. 

Besides, the convective core just appearing in both compo- 
nents is too small (less than one tenth of stellar radius) to in- 
volve a modification of the tracks when changing the overshoot 
parameter a ov from to 0.2. 

Note also that switching from a model atmosphere built as- 
suming an Edding ton grey law or bas ed on Kurucz model at- 
mospheres jKuruc J 1 9791 1 1 992L 1 1 99 j) has no effect as the en- 
velopes of both RS Cha components are fully radiative. 
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4.1 .2. Equation of state: 

Coulomb effects and departure from a perfect gas are very 
small for a 2M star. Hence no variation of the tracks, what- 
ever the evolutionary st age, is observed wh en using either the 
equation of state EFF feggleton et alJll973l) or the equation of 
state OPAL jRogers etalll996h . 

4.1 .3. Burning of light species: 

The light species 2 H, 7 Li and 7 Be burn at the beginning of 
the PMS phase during the first million years of the star, dur- 
ing which the gravitational energy is dominant. The energy re- 
leased by the burning of light species is too small to perturb the 
subsequent evolutionary tracks. 

4.1 .4. Mass loss: 

PMS stars of intermediate mass are known to possess stellar 
winds with mass l oss ra te in the range 1CT 8 - 10~ 7 M .yr~ l (e.g. 
Bouret & Catala They can also accrete mass from a 

surrounding disk with mass accretion rates similar to those of 
the wind. 

We have computed PMS evolutionary tracks assuming 
mass loss or accretion with a rate of 10 _8 M G .yr _1 . Initial 
masses have been chosen such that masses of the models at 
the age of RS Cha (tuscha) are equal to the observed masses 
of RS Cha. As a consequence, although the evolutionary tracks 
differ before and after the age of RS Cha depending on the as- 
sumed mass loss or accretion, the error boxes of RS Cha are 
unaffected. 

4.2. Unchanged luminosity ratio 

We now discuss the physical parameters that affect the position 
of the error boxes in the HR diagramme but do not allow to 
solve the problem of the luminosity ratio. From now on and for 
the sake of visibility, we have removed the evolutionary tracks 
from the figures. Only error boxes obtained from two different 
models along with crosses are kept into the HR diagram. Boxes 
represented by dashed lines are calculated using standard mod- 
els as defined in Sect.2.1. 

4.2.1. Opacity 

We change the opacity (k) in the whole star by diminishing it 
by 10%. Fig.[3]shows that the error boxes are shifted towards 
larger luminosities and temperatures. A decrease in opacity 
makes the energy release easier and in order to maintain the 
pressure balance models become more centrally concentrated 
and hotter which leads to an incre ase in luminosity and effec- 
tive temperature ( Eddinaton Il92dl . However both boxes move 
together towards larger temperatures and luminosities leaving 
the ratio L P /L S unchanged. 
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Fig. 3. Error boxes calculated by decreasing OPAL opacity 
of 10% (full line) and standard OPAL opacity (dashed lines). 
Crosses represent observational error bars in effective temper- 
ature and luminosity (Table for the primary (P) component 
(left) and for the secondary (S) component (right) of RS Cha. 
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Fig. 4. Error boxes calculated by increasing the helium mass 
fraction while keeping Z/X constant: Y = 0.267 (dashed line) 
and Y = 0.300 (full line). 



4.2.2. Helium mass fraction and metallicity 

In paper I, we have determined the metallicity [Fe/H] of RS 
Cha (Table[0, assuming that both stars have the same chemical 
composition and that the metal abundance ratios are the solar 
ones. The knowledge of [Fe/H] gives us a constraint on ~ for 
RS Cha. In order to know fully the composition of the star, we 
have to assume the helium mass fraction Y which is unknown 
for this system. 

We have varied Y between the solar value: Y - 0.267 and 
Y = 0.300 which is a reasonable range compared to the ob- 
served values in various stars and clusters. Fig. 0] shows that 
the luminosity and the effective temperature increase with in- 
creasing Y. An increase of Y leads to an increase of the mean 
molecular weight //. In order to maintain the pressure balance 
the star is more condensed and hotter which results in an in- 
crease in luminosity and effective temperature 
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Fig. 5. Error boxes calculated by decreasing the rate of the 
12 C(p,y) 13 N(e + ,v) 13 C reaction of 10% (full line) from the 
NACRE rate (dashed line) 



Fig. 6. Error boxes calculated by decreasing the abundance of 
carbon: Anc = -0.7 dex (full line), keeping the metallicity con- 
stant. 



We have also varied the metallicity of the system [Fe/H] 
within the error bars. Since heavy elements contribute only for 
a small mass fraction, an increase in heavy element abundances 
does not influence the mean molecular weight fj. very much. 
However an increase in metallicity increases the global opac- 
ity (k). As seen in the previous section, an opacity increase in 
the whole star leads to a decrease in effective temperature and 
luminosity in both stars. 

As both boxes move together in the same direction, we can- 
not reverse the luminosity ratio by varying the helium mass 
fraction or the metallicity. 

4.3. Inverting the luminosity ratio 

So far no parameter, which changes directly or indirectly the 
luminosity transfer into the star, succeeded in inverting the cal- 
culated luminosity ratio of the first component to the second 
one. We now focus on the energy generation. 

The evolutionary stage of RS Cha corresponds to the begin- 
ing of the CNO cycle with the burning of carbon and nitrogen. 
The 12 C— > 14 N chain overheats the internal regions of the star 
leading to the apparition of a convective core responsible for 
the luminosity decrease. The onset of the convective core for 
the primary must be delayed in order to invert the luminosity 
ratio; this means a decrease of the efficiency of overheating of 
the internal regions. This efficiency depends on both nuclear 
reaction rates and the abundances of the species involved in the 
CNO cycle. 

4.3.1. Nuclear reaction rates 

We first decrease the nuclear reaction r ates by 10%, whic h cor- 
responds to their estimated accuracy dAnyulo et al.ll99^ . Fig. 
[5] shows the shift of the boxes resulting from the decrease of 
the rate of the 12 C(p,y) 13 N(e + ,v) 13 C reaction. The luminosity 
of the primary comes closer to that of the secondary, but not 
sufficiently to invert the luminosity ratio. The decrease of the 



rates of the other reactions of the chain 13 C(p,y) 14 N shows a 
similar, although smaller, shift. 

4.3.2. Abundances of Carbon, Nitrogen and Oxygen 

We then decrease the abundance of carbon, nitrogen and oxy- 
gen, changing only the abundance ratios of heavy elements 
while keeping the metallicity unchanged. At the evolutionary 
stage of RS Cha the variation of the abundance of I7 is very 
weak, compared to those of C and N, and the abundance of 
ls O does not vary. The reason is the low rate of the reactions 
16 O(p,y) 17 F0S + v) 17 O and 17 0(p,or) 14 N (NACRE, Angulo et al. 
1999). Hence we do not consider the oxygen abundance any 
further. 

Hereafter we call n, the logarithm of the abundance by 
number TV, of a chemical element i in a scale where the abun- 
dance of hydrogen by number is Nh - 10 12 , that is 

n,- = log ^ + 12 (1) 

Fig. [6] shows the shift of the error boxes resulting from a 
decrease of the carbon abundance nc by 0.7 dex. As we re- 
move some carbon, the reactions I2 C— » 13 C and I3 C— > I4 N be- 
come less efficient. The nuclear energy released e is propor- 
tional to L/m, L and m being the local luminosity and mass; 
the radiative gradient V ' md is proportional to where k, 

P and T are the local o pacity, pressure and temperature (e.g. 
( Kirjrjenhahn & Weiaert 1990)). We find that decreasing the 
abundance of carbon leads to a negligible increase of k, P and 
T in the center of the star, compared to the large decrease of 
the L/m ratio. Therefore, by reducing the carbon abundance, 
the drop of the produced nuclear energy leads to a dominant 
decrease of the L/m ratio, and therefore to a decrease of the 
radiative gradient V ra£ / which delays the onset of the convec- 
tive core. Fig.0shows the variation with time of the radius of 
each star in model 1 (normal carbon abundance) and model 2 
(n c depleted by 0.7 dex). While at ages lower than ~9 Myr, the 
radii of the models as well as the limit between the radiative 
core and the convective envelope are the same, at ages higher 
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than ~9 Myr -after the onset of the CNO cycle- the convective 
core appears later and remains smaller in model 2. 

As a result, the expansion of the internal shells is delayed, 
as well as the decrease in total luminosity (see Sec. 2). The 
evolutionary track is extended towards larger luminosities and 
temperatures and we therefore obtain a luminosity ratio (jj) 
greater than one, as observed. We notice the same phenom- 
ena with the decrease of the nitrogen abundance, which makes 
the reactions 14 N(/?,y) 15 0(/? + , v) 15 N and 15 N(/?,y) 16 less effi- 
cient. 

5. Calibration of RS Cha 

We have finally found that the luminosity ratio is sensitive 
to the abundances of carbon and nitrogen. We now search a 
model matching properly the observations. This model is not 
very different from our standard model. We must change only 
the abundance ratios of heavy species and the helium mass frac- 
tion. Our final models are obtained by decreasing the carbon 
and nitrogen abundances by 0.6 and 0.5 dex respectively. The 
decrease of 0.6 dex and 0.5 dex of the number of atoms of car- 
bon and nitrogen gives new abundances equal to 8 .25 and 7.77, 
respectively (the solar values are 8.55 and 7.971. baflon etalJ 
(2004) measured the abundances of various species in 69 young 
OB stars in 25 clusters. They found abundances comprised be- 
tween 8.21 and 8.49 for carbon and between 7.19 and 7.83 for 
nitrogen, showing that our values agree well with abundances 
observed in young stars. We should also note that our results 
are going in the same direction as the recent solar abundances 
determination of Asplund et al. (2004,2005a,b). They found 
a decrease of 0.16 dex of the carbon abundance and of 0.11 
dex of the nitrogen abundance and of 0.21 dex of the oxy gen 
abundance, compared to the abundances of iGrevesse & Noelsl 
Jl993h used in the CESAM code. 

Having the appropriate luminosity ratio between both com- 
ponents of the system, we need to increase the helium mass 
fraction in order to match the calculated parameters of RS Cha 
to the observed ones (Table [Q. The error bars in masses and 
radius of the models are plotted in Fig. [8] We find a value 
of Y = 0.272 (log(N He /N H ) + 12 = 10.99), reasonable com- 
pared to the helium mass fraction observed in young stars 
(iMathvs et all2002l) . 

Finally, it remains to be checked that both stars have the 
same age. We find that a common interval of ages does indeed 
exist inside the error boxes indicating that our final models for 
both components can have the same age. We illustrate this by 
plotting the 9.50 Myr isochrone crossing both boxes in Fig. 
00 T his age is of the same order as previous d eterminations 
(Mamaie k et al.ll2000t iLuhman & Steeghsl2004l) . However we 
stress our age determination is based on an observed metallicity 
contrary to previous works. 

Besides, in a binary system the ratio of the effective temper- 
atures is better constrained than the abso lute values of the tem- 
peratures of the individual stars. Idausen & Nordstrom1(ll980l) 
measured very accurately the photometric value of this ratio : 
Tp/Ts = 1.055 + 0.004. In order to take into account this ad- 
ditionnal constraint, we checked that our calibration model can 



reproduce simultaneously this ratio and the observed radii of 
both stars, by plotting error boxes in the radius versus Ip/Ts 
diagram as follows. We first calculated the effective tempera- 
tures ranges delimited by the error boxes in masses and radii of 
FigODfor the primary and secondary components on the same 
isochrone. Then we restricted these intervals to the primary 
and secondary temperatures which satisfy the observed ratio 
of Tp/Ts within the error bars. Finally we calculated by lin- 
ear interpolation the corresponding radii ranges. We obtain two 
boxes as shown in the lower right panel of the Fig. [8] We find 
that our calibration models reproduce the observed radii of both 
stars as well as the observed effective temperature ratio of the 
system. 

We have seen that the small area of the PMS tracks where 
RS Cha is observed is only weakly sensitive to various physical 
inputs and corresponding stellar models are then only weakly 
constrained. Fortunately more severe constraints can be ob- 
tained by a seismic sounding of the stellar interior, provided 
the stars oscillate. 

6. RS Cha: a binary system of pulsating PMS stars 

6.1. Observations of delta-scuti type pulsations 

Many authors jAndersenl Il975t iMdnallv & Austin! Il977t 
IPalla & StahleJ200ll) discussed the possibility of the presence 
of oscillations in one of the components of RS Cha but re- 
ported no direct observations of these oscillations. Andersen 
(1975) mentioned hints of variability in the residuals from the 
observed radial velocity curve of the primary but could not be 
more conclusive with the data at his disposal. Marconi & Palla 
(1998) calculated the theoretical location of the instability strip 
of PMS stars in the HR diagram for radial modes. Palla & 
Stahler (2001) located RS Cha in the HR diagram and found 
that the secondary component lies inside the instability strip 
and ought to be a PMS 6 Scuti star. The primary component 
was found outside but near the left of the instability strip for 
fundamental radial mode which cannot exclude that the star is 
pulsating in radial overtones and/or nonradial modes. 

Finally, using our data, we have recently shown, that both 
components of RS Cha are probably pulsating (paperl). We 
have pointed out temporal variations in the residuals from ra- 
dial velocity curves of both stars. These variations may be peri- 
odic with periods around one hour. We therefore ascribed them 
to 6 Scuti type oscillations. 

Unfortunately our data were not accurate enough and did 
not cover sufficient time to determine precisely the periods. 
However we can wonder whether the period around one hour 
found with our data belongs to the period range of the excited 
modes that can be expected theoretically. 

6.2. Theoretical oscillation modes and periods 

Theoretical periods of pulsation modes of both compo- 
nents of RS Cha w ere calculated using Dupret's code MAD 
JPupret et al.ll2005T) . This nonadiabatic code includes a time- 
dependent interaction between pulsation and turbulent convec- 
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Fig. 7. Radius of the primary (left) and secondary (right) star plotted as a function of the age of the standard model 1 in thick- 
dashed-line, and model 2 built with a decrease of 0.7 dex of carbon in thick-full-line. Before ~ 10 Myr both radii are superimposed. 
The thin lines represent the limit between radiative and convective zones in model 1 (thin-dashed-line) and in model 2 (thin-full- 
line). The two vertical full lines represent the extrema of the age of the primary and secondary components of RS Cha (RS Cha 
P and RS Cha S) calculated from the error boxes P and S of Fig. [6] We see that the convective core in the secondary component 
is less developped than that of the primary one 



tion dGrigahcene et alJl2005l) . The blue and red edges of the 
instability strip are accurately determined with thi s code for ra- 
dial and non-radial modes of MS and PMS stars JPupret et alJ 
EooilCrigahcene et all2006h. 

The stellar parameters of the selected models for RS Cha 
components are those resulting from our calibration (see Sec. 
5). 

In paper I, it was shown that the system is synchronized and 
the orbital period of the system has been determined accurately, 
we therefore know the rotational period of both stars: P rot = 
1 .67 d. For sake of simplicity we assume a uniform rotation 
with the surface angular velocity Q. = 2nlP rot . This enables us 
to compute the rotational splitting for each mode , according to 
the first order perturbation theory dLedouxll9 58i 

Vnbn = mO(l - C nt ) (2) 

in the observer reference frame. C n ,t is the Ledoux constant. 
These m + components of the multiplets are included in the 
frequency set and they are assumed stable or unstable accord- 
ing to the same status of the associated m = mode. 

Fig- shows the frequency interval over which modes are 
unstable and therefore expected to be detected. Comparison be- 
tween the plots of each component illustrates a well known 
result for 5 Scuti like pulsations: unstable modes have higher 
frequencies and radial order towards the blue side of the insta- 
bility strip (here the primary component). 

We find that the periods of excited modes are around one 
hour for both stars (those of the secondary are a little longer 
than in the primary). On one hand, this supports the conclusion 
that the observed period corresponds to a low order pressure 
mode of delta-Scuti type and on the other hand this shows our 
models are not too far from reality. 

In Sect. 4.1.1., we have seen that including convective core 
overshoot ot not does not change the location of the RS Cha 



models in the HR diagram indicating that the structures of the 
models are not too much affected when considering these pho- 
tometric observables. We now consider the effet of including 
overshoot on the frequencies of each component of the system 
separately as an indication of the efficiency of the seismic di- 
agnostic ability of RS Cha system. 

For each component of RS Cha, we compare the frequen- 
cies for each mode nl computed from two models which dif- 
fer by the value of the overshoot parameter: one model has no 
overshoot d ov = and the other one is computed assuming 
d ov = Q.2H P . These models selected for comparison have the 
same mass and are located quite close to each other in the HR 
diagram. Fig^|shows frequency differences 



as a function of v / . Frequencies v V refer to frequencies for 
each oscillation mode (n, I) computed from the stellar model 
with convective core overshoot and v , those for the stellar 
model computed assuming no convective core overshoot. The 

ratio -^t; is a scaling factor which corrects for the fact that 

V 1.0 

the models with and without overshoot do not have the same 
radius. 

Differences at high frequency for the secondary arise be- 
cause the selected models with and without overshoot for this 
component have not quite the same radius and this is not com- 
pletely corrected by the scaling v^/v^. . On the other hand 
at low frequency, differences are large for the primary as the 
models which have been chosen for comparison have the same 
radius but not quite the same evolutionary status. FigllOlindi- 
cates that large frequency differences arise for modes which 
are in the unstable range for both stars and if detected could be 
quite discriminant. 
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Fig. 8. Evolutionary tracks (dash-dotted lines) and error boxes (full lines) calculated by decreasing the abundances of carbon and 
nitrogen: Anc = -0.6 dex and An^ = -0.5 dex. The thick line is the 9.50 Myr isochrone. Crosses represent observational error 
bars in effective temperature and luminosity (Table 1) for the primary (P) component (left) and for the secondary (S) component 
(right) of RS Cha. The lower right pan el shows the constra ints on the radii of both stars when the photospheric measure of the 
effective temperature ratio of IClausen & Nordstrom! (Tl980) is taken into account. Crosses represent the error bars in radii and 
temperature ratio for both stars. Each line in the hatched area represents the allowed values of R for a given value of Tp/Ts- 
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Fig. 9. Theoretical frequencies of pulsation modes obtained with non-adiabatic calculations, for the primary (left) and the secondary (right) 
stars. Each bar represents one mode as a function of its frequency (in fiHz) and its degree I (vertical direction). Dashed bars represent stable 
modes and solid bars unstable modes. Label p is for p mode and label g for g mode. Numbering of the label p and g referes to the radial order 
of the mode 
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Fig. 10. Frequency differences as a function of the frequency for each 
oscillation mode (n, £) between frequencies of one stellar model with 
overshoot and one without overshoot for the primary RS Cha compo- 
nent on the right and the secondary RS Cha component on the left. 
Black (resp.open) dots represent unstable (resp. stable) modes 



7. Conclusions 

In this paper, we carried out a stellar modelling of a binary 
system whose physical parameters - masses, radii, luminosi- 
ties and effective temperatures are well constrained by obser- 
vations. The RS Cha system is composed of two similar inter- 
mediate mass stars in a PMS stage of evolution quite close to 
the ZAMS at the onset of CNO burning. Assuming that both 
stars have the same chemical composition, a first modelling us- 
ing the spectroscopically determined iron abundance for this 
system and assuming solar abundance ratios could not satisfy 
the observational constraints simultaneously for both stars. 

In order to obtain a consistent modelling of both stars, we 
investigate the sensitivity of the models in the evolutionary 
stage of the RS Cha system - at the end of the PMS phase - 
to various uncertainties in the modelling, including uncertain- 
ties in the input physics and initial abundances. We found that 
the stellar structure is only weakly sensitive to uncertainties in 
the stellar physical description. 

The calibrated models i.e. the final ones which satisfy all 
the observational constraints were obtained only by changing 
the abundance ratios of carbon and nitrogen which then are in 
agreement with known abundances of other young stars. Once 
the chemical composition is set as well as the other stellar pa- 
rameters intrinsic to RS Cha stars, the error bars on the ef- 
fective temperatures and luminosities, although already quite 
small, remain still too large to put the standard physical input 
of our stellar models at default. These models also satisfy the 
constraints on the effective temperature ratio. 

We find that there is no need to include rotation in our mod- 
els as they are able to reproduce the observed properties of RS 
Cha within their error bars. We note that similar conclusions 



are drawn in their study of EK Cep binary sys- 

tem. Furthermore we have checked that the effects of rotation 
on the global parameters in the HR dia gram (L, r P ff) estimated 
from the grid of rotating atmospheres of Perez Hernandez et al. 
dl999l) for the vsin/ values of the components (Table 1) are 
negligible with respect to the error bars. 

In order to obtain useful constraints on the physical de- 
scription in this evolutionary phase, we therefore need to find 
other properties of stars than the photometric and spectroscopic 
ones, that-is seismic properties. The binary system RS Cha is 
a good example as we have recently proved that both compo- 
nents are pulsating. The detected variations indicate periodic- 
ities around 1 hour which we find to fall within the range of 
theoretically opacity driven modes with radial orders P3-P4 for 
the primary and P4-P5 for the secondary. We find that in this 
range the frequency spectra for both stars are regular when in- 
cluding modes with spherical degrees I — 0, 2, an advantage 
when dea ling with mode id entification and possible seismic in- 
ferences dSuran e t al. 2001). 

For a deeper study however, more observations are needed 
in order to determine the seismic properties of both stars more 
precisely. This is not an easy task. This binary system is not 
resolved even with our biggest telescopes. We cannot there- 
fore obtain separate photometric light curves. We can only ob- 
serve the spectroscopic variations of the radial velocities of 
both stars. Getting modes and periods by this way, requires 
long-term multi-site monotoring (which is under preparation). 
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